Generalised fft-ifft structure based frequency division multiplexing transceiver

ABSTRACT

A generalized frequency division multiplexing (GFDM) transceiver system includes a low complex GFDM transmitter with multiple sub-carriers and timeslots having IFFT based modulator for modulating data corresponding to a particular timeslot and different sub-carriers to corresponding sub-carrier frequencies and thereby generating transmittable GFDM data signal, multipath frequency selective fading GFDM channel having uncorrelated channel coefficients corresponding to different paths for transmitting the modulated GFDM data signal and a low complex GFDM receiver configured to operate with said multipath frequency selective fading channel involving channel equalization followed by self-interference equalization to receive the transmitted modulated GFDM data signal and thereby de-modulate the GFDM data signal to obtain the data.

CROSS-REFERENCE TO RELATED APPLICATION(S)

This application claims the benefit under 35 USC § 119(a) of Indian Patent Application No. 2017-31045300, filed on Dec. 16, 2017, in the Indian Intellectual Property Office, the entire disclosure of which is incorporated herein by reference for all purposes.

BACKGROUND Field of the Invention

The present invention relates to Generalized Frequency Division Multiplexing (GFDM). More specifically, the present invention is directed to develop a GFDM transceiver with least computational complexity facilitating implementation of cheap and fast GFDM based communication infrastructure. The present GFDM transceiver can strengthen 5G waveform candidature of the GFDM as it reduces cost of the GFDM based modem and also increases the processing speed (which reduces overall latency).

Background of the Invention

The fifth generation (5G) communication systems aim to cater to a wide range of application with varied requirements [Ref J. G. Andrews et al., “What Will 5G Be?,” IEEE Journal on Selected Areas in Communications 32, no. 6 (June 2014): 1065-82]. Orthogonal Frequency Division Multiplexing (OFDM), has been the celebrated waveform for fourth generation (4G) cellular systems due to its low complexity implementation and frequency selective channel combating feature. However, it is shown to fall behind in terms of 5G waveform requirements, such as, very low out of band (OoB) emission, low latency, immunity to carrier frequency offset (CFO) [Ref: P. Banelli et al., “Modulation Formats and Waveforms for 5G Networks: Who Will Be the Heir of OFDM?: An Overview of Alternative Modulation Schemes for Improved Spectral Efficiency,” IEEE Signal Processing Magazine 31, no. 6 (November 2014): 80-93]. In recent years, many new waveforms have been suggested for 5G and GFDM [Ref N. Michailow et al., “Generalized Frequency Division Multiplexing for 5th Generation Cellular Networks,” IEEE Transactions on Communications 62, no. 9 (September 2014): 3045-61] is one the main contenders among many candidate waveforms for 5G.

GFDM is a block-based waveform which encompasses multiple time and frequency slots. Circular pulse shaping is used to restrict the signals length within its duration and enable the use of cyclic prefix (CP) for combating frequency selective wireless channel. Furthermore, GFDM is shown to have better OoB characteristics than OFDM and can achieve even excellent OoB characteristics with the use of time domain windowing without any spectral efficiency loss, which qualifies it to be a contending waveform for cognitive radio applications.

The GFDM receivers in multipath channel can be broadly categorized as (a) two-stage receiver [Ref: N. Michailow et al., “Bit Error Rate Performance of Generalized Frequency Division Multiplexing,” in 2012 IEEE Vehicular Technology Conference (VTC Fall), 2012, 1-5] and (b) one-stage receiver [Ref: Matthe et al., “Reduced Complexity Calculation of LMMSE Filter Coefficients for GFDM.”]. In a two-stage receiver, channel equalization is followed by GFDM demodulation, while in the one stage receiver, the effect of channel and GFDM modulation is jointly equalized. In the two stage GFDM receiver, channel equalization can be implemented using low computational load. The second stage, which is a self-interference equalizer, can be implemented using linear or non-linear receiver.

If M and N represent number of time and frequency slots respectively in GFDM based communication, the implementation of the transmitter, Matched Filter (MF) receiver (self-interference equalizer) and Zero-Forcing (ZF) receiver (self-interference equalizer) involves a complexity of O(M²N²) [Ref: Shashank Tiwari, Suvra Sekhar Das, and Kalyan Kumar Bandyopadhyay, “Precoded Generalised Frequency Division Multiplexing System to Combat Inter-Carrier Interference: Performance Analysis,” IET Communications, Sep. 10, 2015,] while the complexity of Minimum Mean Square Error (MMSE) receiver (self-interference equalizer) is O(M³N³). When N˜10³s and M˜10s, the count of computations becomes very high. This high complexity hinders practical implementation of GFDM transceivers. It is known that the Joint-MMSE receiver outperforms two-stage receivers [Ref: Michailow et al., “Generalized Frequency Division Multiplexing for 5th Generation Cellular Networks.]. Despite good BER properties, Joint-MMSE receiver is the most complex linear receiver for constituting GFDM receiver since it involves large matrix multiplications, inversion and O(M³N³) computations.

Some attempts to reduce the complexity of GFDM transmitter and two stage receiver is reported in the past. The sparsity of prototype pulse shape in frequency domain is exploited to design a low complexity transmitter in N. Michailow et al [Ref: N. Michailow et al., “Generalized Frequency Division Multiplexing: Analysis of an Alternative Multi-Carrier Technique for next Generation Cellular Systems,” in 2012 International Symposium on Wireless Communication Systems (ISWCS), 2012, 171-75,] and a low-complexity MF receiver in I. Gaspar et al. [Ref: I. Gaspar et al., “Low Complexity GFDM Receiver Based on Sparse Frequency Domain Processing,” in Vehicular Technology Conference (VTC Spring), 2013 IEEE 77th, 2013, 1-6,]. The complexity is reduced to O(MN log₂ (MN)+MN²) but it comes with increase in BER. Behrouz and Hussein proposed frequency spreading based GFDM transmitter in Behrouz Farhang-Boroujeny et al [Ref: Behrouz Farhang-Boroujeny and Hussein Moradi, “Derivation of GFDM Based on OFDM Principles,” in 2015 IEEE International Conference on Communications (ICC), 2015, 2680-85] based on the principles of frequency spreading filter bank multi carrier (FMBC) transmitter proposed in M. Bellanger [Ref: M. Bellanger, “Physical Layer for Future Broadband Radio Systems,” in 2010 IEEE Radio and Wireless Symposium (RWS), 2010, 436-39,]. The complexity of the transmitter is O(MN log₂(N)+M²N). Periodicity of complex exponential is exploited in Hao Lin et al [Ref: Hao Lin and Pierre Siohan, “Orthogonality Improved GFDM with Low Complexity Implementation,” in 2015 IEEE Wireless Communications and Networking Conference (WCNC), 2015, 597-602,] and Maximilian Matthé et al. [Ref: Maximilian Matthé et al., “Precoded GFDM Transceiver with Low Complexity Time Domain Processing,” EURASIP Journal on Wireless Communications and Networking 2016, no. 1 (May 25, 2016): 1] to attain O(MN log₂(N)+M²N) complexity of GFDM transceivers. Similar order of complexity is achived by using block circulant property of multiplication of modulation matrix and its Hermitian in A. Farhang et al [Ref: A. Farhang, N. Marchetti, and L. E. Doyle, “Low-Complexity Modem Design for GFDM,” IEEE Transactions on Signal Processing 64, no. 6 (March 2016): 1507-18].

Some attempts to reduce the complexity of Joint-MMSE receiver also have been reported in the recent past, there is either limited gain achieved or complexity gain is traded off with BER loss. Matthe et. al. [Ref: Matthe et al., “Reduced Complexity Calculation of LMMSE Filter Coefficients for GFDM] exploited the block circulant property of the modulation matrix to achieve the complexity of O(MN² log M). Authors in M. Matthe, D. Zhang, and G. Fettweis,” Iterative Detection Using MMSE-PIC Demapping for MIMO-GFDM Systems,” in European Wireless 2016; 22th European Wireless Conference, 2016, 1-7 have exploited the sparsity of prototype pulse in the frequency domain to achieve the complexity of O(M³N). Recently, authors in Zhang et al., “A Study on the Link Level Performance of Advanced Multicarrier Waveforms Under MIMO Wireless Communication Channels,” IEEE Transactions on Wireless Communications 16, no. 4 (April 2017): 2350-65, have computed MMSE equalization of FFT of data vectors while exploiting the sparsity of prototype pulse in the frequency domain to achieve the complexity of O(MN log₂MN).

It is thus there has been a need for developing a new simple GFDM transceiver by reducing the complexity of the GFDM transceiver without any significant loss in the BER performance.

It is thus basic object of the present invention is to develop a GFDM transceiver which would have the least computational complexity when compared with other GFDM transceiver structures without exhibiting any significant loss in the BER performance.

Another object of the present invention is to develop a GFDM transceiver which would enable implementation of GFDM based communication cheaper in terms of cost and faster in terms of processing speed.

Another object of the present invention is to develop a GFDM transceiver which would include digital signal processing blocks such as Fast Fourier Transform (FFT), Multiplier, Adder and like having a novel interaction there between for providing a novel signal propagation path from one end to the other end.

Yet another object of the present invention is to develop a GFDM transceiver which would be adapted to implement GFDM based communication in fifth generation cellular systems, machine-type communication, Internet of Things (IOT), Tactile Internet, Cognitive Radio etc.

SUMMARY

Thus, according to the basic aspect of the present invention there is provided a generalized frequency division multiplexing (GFDM) transceiver system comprising

low complex GFDM transmitter with multiple sub-carriers and timeslots having IFFT based modulator for modulating data corresponding to a particular timeslot and different sub-carriers to corresponding sub-carrier frequencies and thereby generating transmittable GFDM data signal;

multipath frequency selective fading GFDM channel having uncorrelated channel coefficients corresponding to different paths for transmitting the modulated GFDM data signal; and

low complex GFDM receiver configured to operate with said multipath frequency selective fading channel involving channel equalization followed by self-interference equalization to receive the transmitted modulated GFDM data signal and thereby de-modulate the GFDM data signal to obtain the data.

In a preferred embodiment of the present generalized frequency division multiplexing (GFDM) transceiver system, the low complex GFDM transmitter includes

N-point IFFT operator to receive the data corresponding to a particular timeslot and different sub-carriers and modulate the same to corresponding sub-carrier frequencies;

means for shuffling physical connections in the N-point IFFT operator for grouping the modulated data to sub-carrier numbers, whereby, in each group, the data is converted into frequency domain using M-point FFT operator and multiplied with a precomputed weight and thereafter converted back into time domain by using M-point IFFT operator;

means for shuffling physical connections in the M-point IFFT operator for grouping the data according to time slots and generate transmittable GFDM data signal.

In a preferred embodiment of the present generalized frequency division multiplexing (GFDM) transceiver system, the low complex GFDM transmitter generated transmittable GFDM data signal is

${{x\lbrack n\rbrack} = {\frac{1}{\sqrt{N}}{\sum\limits_{m = 0}^{M - 1}\; {\sum\limits_{k = 0}^{N - 1}{d_{m,k}{g\left\lbrack {n - {mN}} \right\rbrack}_{MN}\mspace{11mu} e^{\frac{j\; 2\; \pi \; {nk}}{N}}}}}}}\;$

for N sub-carriers and M timeslots, where g(n), n=0, 1, . . . , MN−1 is MN length filter response and d_(m,k) ∈C, m=0, 1, . . . , M−1, k=0, 1, . . . , N−1 is QAM modulated data symbol.

In a preferred embodiment of the present generalized frequency division multiplexing (GFDM) transceiver system, the GFDM transmitted signal is critically sampled Inverse Discrete Gabor Transform (IDGT) of d by using the IDGT matrix factorization whereby Modulation Matrix, A can be given as,

$A = {{G \times U_{N}} = {\begin{pmatrix} \Psi_{0} & \Psi_{M - 1} & \ldots & \Psi_{1} \\ \Psi_{1} & \Psi_{0} & \ldots & \Psi_{2} \\ \vdots & \vdots & \ddots & \vdots \\ \Psi_{M - 1} & \Psi_{M - 2} & \ldots & \Psi_{0} \end{pmatrix}\begin{pmatrix} W_{N} & \; & \; \\ \; & \ddots & \; \\ \; & \; & W_{N} \end{pmatrix}}}$

where, Ψ_(m)=diag{g[mN], g[mN+1], . . . , g[mN+N−1]} for 0≤m≤M−1, is N×N diagonal matrix and W_(N) is N×N normalized IDFT matrix

In a preferred embodiment of the present generalized frequency division multiplexing (GFDM) transceiver system, the uncorrelated channel coefficients corresponding to different paths for transmitting the modulated GFDM data signal constitutes channel impulse response vector given as h=[h₀, h₁, . . . h_(L-1)]^(T) where L is channel length and h_(i), for 0≤i≤L−1, represents complex baseband channel coefficient of (i+1)^(th) path, which is assumed to be zero mean circular symmetric complex Gaussian whereby received vector of length N_(CP)+NM+L−1 (for N_(cp)≥L) is given by,

Z _(cp) =h*X _(cp) +v _(cp)

where V_(cp) is AWGN vector of length MN+N_(CP)+L−1 with elemental variance σ² _(v).

In a preferred embodiment of the present generalized frequency division multiplexing (GFDM) transceiver system, the low complex GFDM receiver includes two staged receiver or Joint-MMSE Receiver whereby, the data obtained from the received GFDM data signal by involving channel equalization followed by self-interference equalization.

In a preferred embodiment of the present generalized frequency division multiplexing (GFDM) transceiver system, the two-staged receiver includes two staged receiver includes:

M-point FFT operator for grouping channel equalized received GFDM data signal according to sub-carrier numbers followed by

-   -   converting samples of each group into frequency domain using the         M-point FFT operator;     -   multiplying the converted samples with pre-computed weights by         multiplier means;     -   converting back the multiplied samples into time domain by using         M-point IFFT operator;

means for shuffling physical connections in the M-point IFFT operator for regrouping the converted samples according to time slots followed by converting samples of each group into frequency domain using N-point FFT operator;

multiplying the converted samples with pre-computed weights by multiplier means to obtain self-interference equalized data signal.

In a preferred embodiment of the present generalized frequency division multiplexing (GFDM) transceiver system, the joint-MMSE receiver includes

MN-point FFT operator to convert the equalized received GFDM data signal into frequency domain;

means to multiply the received signal in frequency domain with complex valued channel information-based weights and then converting back to time domain using MN-point IFFT operator;

means for reshuffling the physical connections in the MN-point IFFT operator for grouping the time domain converted samples according to subcarrier number followed by

-   -   converting the samples of each group into frequency domain using         the M-point FFT operator;     -   multiplying the converted samples with pre-computed weights;

means for reshuffling the physical connections in the M-point FFT operator for regrouping the converted samples according to time slots followed by

converting samples of each group into frequency domain using the N-point FFT operator;

processing the converted samples following Algorithm 1;

means for reshuffling the physical connections in the N-point FFT operator for regrouping the processed samples according to sub-carrier number followed by

-   -   converting the samples of each group to time domain using the         M-point IFFT operator;     -   multiplying the converted samples with complex weights computed         using Algorithm 2;

means for reshuffling the physical connections in the M-point IFFT operator for regrouping the multiplied samples according to time slots to obtain equalized MN-point samples.

In a preferred embodiment of the present generalized frequency division multiplexing (GFDM) transceiver system, the Algorithm 1 enabling development of the low complexity receiver structure including low complexity multiplication to obtain ρ=Eκ using Taylor Series expansion.

In a preferred embodiment of the present generalized frequency division multiplexing (GFDM) transceiver system, the Algorithm 2 enabling development of the low complexity receiver structure including low complexity multiplication to obtain ρ=Eκ using CG method.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows a wireless communication system involving GFDM based transceiver in accordance with an embodiment of the present invention.

FIG. 2 shows baseband block diagram of GFDM Transceiver in multipath fading channel with two staged receiver in accordance with an embodiment of the present invention.

FIG. 3 shows low complexity implementation of GFDM transmitter in accordance with an embodiment of the present invention.

FIG. 4 shows low complexity implementation of GFDM receiver in multipath fading channel in accordance with an embodiment of the present invention.

FIGS. 5A to 5D show ratio of power in constituent matrices of Φ_(u) i.e.

$\frac{{\Delta_{u}}_{F}}{{Z_{u}}_{F}}$

averaged over h and u∈[0 M−1] i.e. ρ (in dB). (Raised Cosine (RC) pulse shape is considered).

FIG. 6 shows implementation of present low-complexity joint-MMSE GFDM receiver in accordance with an embodiment of the present invention.

FIG. 7 shows BER Performance of proposed GFDM transceiver and direct implementation GFDM two stage transceiver for 16 QAM modulation in fading channel. N=128, M=8.

FIGS. 8A and 8B show BER Performance of proposed GFDM transceiver for (a) Roll off factor 0.3 and (b) Roll off factor 0.9.

DETAILED DESCRIPTION

In this specification, following notations are used. Vectors are represented by bold small letters (x), matrices are represented by bold capital letters (X) and scalars are represented as normal small letters (x). I_(N) represents identity matrix with order N and j=√{square root over (−1)}. W_(L) represents L-order normalized IDFT matrix. Kronecker product operator is given by ⊗. diag{.} is a diagonal matrix whose diagonal elements are formed by the elements of the vector inside or diagonal elements of the matrix inside. circ{.} is a circulant matrix whose first column is given by the vector inside. E_(h){.} is expectation of expression inside with respect to random vector h. The round-down operator *. , rounds the value inside to the nearest integer towards minus infinity. The superscripts (.)^(T) and (.)^(H) indicate transpose and conjugate transpose operations, respectively.|.| operator computes absolute value of elements inside. trace{.} computes the trace of matrix inside. ∥.∥ is Frobenous norm of matrix inside. FFT_((.)) and IFFT_((.)) denote (.)-point FFT and IFFT respectively.

GFDM Transceiver System:

The accompanying FIG. 1 shows a block representation of the wireless communication system involving GFDM based transceiver in accordance with an embodiment of the present invention. As shown in the referred figure, the present transceiver system includes GFDM transmitter, wireless channel and GFDM receiver.

Transmitter:

The present invention discloses a GFDM system with N sub-carriers and M timeslots. The MN length prototype filter is g(n), n=0, 1, . . . , MN−1. QAM modulated data symbol is d_(m,k)∈C,m=0, 1, . . . , M−1, k=0, 1, . . . , N−1. It is assumed that data symbols are independent and identical i.e. E[d_(m,k)d_(m′,k′)*]=σ_(d) ²δ_(m-m′,k-k′). The transmitted GFDM signal can be written as,

$\begin{matrix} {{x\lbrack n\rbrack} = {\frac{1}{\sqrt{N}}{\sum\limits_{m = 0}^{M - 1}\; {\sum\limits_{k = 0}^{N - 1}\; {d_{m,k}{g\left\lbrack {n - {mN}} \right\rbrack}_{MN}{e^{\frac{j\; 2\; \pi \; {nk}}{N}}.}}}}}} & (1) \end{matrix}$

The transmitted signal can also be written as,

x=A _(MN×MN) d _(MN×1),  (2)

where d=[d₀ d₁ . . . d_(M-1)]^(T) is the data vector, where d_(m)=[d_(m,0) d_(m,1) . . . d_(m,N-1)]^(T), where, m=0, 1 . . . M−1, is the N length data vector for m^(th) time slot and A is the modulation matrix which can be given as,

A=[g M ₁ g . . . M _(N-1) g|T ₁ g T ₁ M ₁ g . . . T ₁ M _(N-1) g| . . . |T _(M-1) M ₁ g . . . T _(M-1) M _(N-1) g],  (3)

where, g=[g[0] g[1] . . . g[MN−1]]^(T) is MN length vector which holds the prototype filter coefficients,

${M_{l}{g\lbrack n\rbrack}} = {{g\lbrack n\rbrack}e^{\frac{j\; 2\; \pi \; \ln}{N}}}$

is the modulation operator and T_(r)=g(n−rN)_(MN) is the cyclic shift operator.

CP of length N_(CP) is prepended to x. After adding CP, transmitted vector, x_(cp), can be given as,

x _(cp)=[x(MN−N _(cp)+1:MN);x].  (4)

Channel:

Let, h=[h₀, h₁, . . . h_(L-1)]^(T) be L length channel impulse response vector, where, h_(i), for 0≤i≤L−1, represents the complex baseband channel coefficient of (i+1)^(th) path [27], which is assumed to be zero mean circular symmetric complex Gaussian (ZMCSC). It is also assumed that channel coefficients related to different paths are uncorrelated. It is considered, N_(cp)≥L. Received vector of length N_(cp)+NM+L−1 is given by,

z _(cp) =h*x _(cp) +v _(cp),  (5)

where v_(cp) is AWGN vector of length MN+N_(cp)+L−1 with elemental variance σ_(v) ².

Receiver:

The first N_(cp) samples and last L−1 samples of y_(cp) are removed at the receiver i.e. y=[y_(cp)(N_(cp)+1: N_(cp)+MN)]. Use of cyclic prefix converts linear channel convolution to circular channel convolution when N_(cp)≥L. The MN length received vector after removal of CP can be written as,

z=HAd+v,  (6)

where H is circulant convolution matrix of size MN×MN and v is WGN vector of length MN with elemental variance σ_(v) ². Since H is a circulant matrix, y can be further written as,

z=W _(MN) ΛW _(MN) ^(H) Ad+v,  (7)

where, A=diag{{tilde over (h)}(0), {tilde over (h)}(1) . . . {tilde over (h)}(MN−1)} is a diagonal channel frequency coefficients matrix whose r^(th) coefficient can be given as,

${\overset{\sim}{h}(r)} = {\sum\limits_{s = 0}^{L - 1}\; {{h(s)}e^{\frac{j\; 2\; {\pi s}\; r}{MN}}}}$

where, r=0,1 . . . MN−1.

In this invention, two stage as well as one stage receiver is considered.

Two Stage Receiver.

For two stage receiver, channel equalized vector, y, can be given as,

$\begin{matrix} {{y = {{W_{MN}\Lambda_{eq}W_{MN}^{H}z} = {{aAd} + b + v}}},} & (8) \\ {{where},\mspace{14mu} {\Lambda_{eq} = \left\{ \begin{matrix} {\; \Lambda^{- 1}} & {{for}\mspace{14mu} {ZF}\mspace{11mu} {FDE}} \\ \left\lbrack {{\Lambda^{H}\Lambda} + {\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I_{MN}}} \right\rbrack^{- 1} & {\Lambda^{H}{for}\mspace{14mu} {MMSE}\mspace{11mu} {FDE}} \end{matrix} \right.}} & \; \\ {{where},\; {v = {W_{MN}\Lambda_{eq}W_{MN}^{H}v}},} & \; \\ {a = \left\{ {{\frac{1}{MN}{\sum\limits_{r = 0}^{{MN} - 1}\; {\frac{\begin{matrix} {{1\mspace{14mu} {for}\mspace{14mu} {ZF}} - {FDE}} \\ {{{\overset{\sim}{h}(r)}^{2}}} \end{matrix}}{{{\overset{\sim}{h}(r)}^{2}{+ \frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}}}}{for}\mspace{14mu} {MMSE}\mspace{11mu} {FDE}}}},} \right.} & \; \end{matrix}$

b is residual interference, given in (9) and v=W_(MN)Λ_(eq)W_(MN) ^(H)v is post-processing noise.

$\begin{matrix} {b = \left\{ \begin{matrix} {0\mspace{14mu} {for}\mspace{14mu} {ZF}\mspace{11mu} {FDE}} \\ {{Ad}\mspace{14mu} {for}\mspace{14mu} {MMSE}\mspace{11mu} {FDE}} \end{matrix} \right.} & (9) \end{matrix}$

Channel equalized vector, y, is further equalized to remove the effect of self-interference. Estimated data, d, can be given as,

d=A _(eq) y,  (10)

where, A_(eq) is GFDM equalization matrix which can be given as,

$A_{eq} = \left\{ \begin{matrix} {A^{H}\mspace{11mu} {for}\mspace{14mu} {MF}\mspace{14mu} {Equalizer}} \\ {A^{- 1}\mspace{14mu} {for}\mspace{14mu} {ZF}\mspace{14mu} {Equalizer}} \\ {\left\lbrack {R_{v} + {A^{H}A}} \right\rbrack^{- 1}A^{H}} \\ {{for}\mspace{14mu} {biased}\mspace{14mu} {MMSE}\mspace{14mu} {Equalizer}} \\ {\;^{\Theta}{{gfdm}^{- 1}\left\lbrack {R_{v} + {A^{H}A}} \right\rbrack}^{- 1}A^{H}} \\ {{{for}\mspace{14mu} {unbiased}\mspace{14mu} {MMSE}\mspace{14mu} {Equalizer}},} \end{matrix} \right.$

for unbiased MMSE Equalizer, where, R_(v)=E[vv^(H)] is noise correlation matrix after channel equalization. In the case of AWGN,

$R_{v} = {\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}{I.}}$

For multipath fading channel, R_(v) is a full matrix since the noise after channel equalization is colored. Θ_(gfdm) ⁻¹ is a diagonal bias correction matrix for GFDM-MMSE equalizer, where,

$\Theta_{gfdm} = {{diag}{\left\{ {\left\lbrack {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I} + {A^{H}A}} \right\rbrack^{- 1}A^{H}A} \right\}.}}$

One-Stage Receiver (Joint-MMSE Receiver):

Joint-MMSE equalizer can be of two types, namely, (1) biased-Joint MMSE and (2) unbiased-Joint-MMSE. Equalized data symbol vector, d_(JP), can be given as, d_(JP)=B_(eq)y, where, B_(eq) is Joint-MMSE equalizer matrix and can be given as,

$\begin{matrix} {B_{eq} = \left\{ \begin{matrix} {{\left\lbrack {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I} + {({HA})^{H}{HA}}} \right\rbrack^{- 1}({HA})^{H}\mspace{14mu} {for}\mspace{14mu} {biased}\mspace{14mu} {Joint}} - {{MMSE}\mspace{14mu} {Equalizer}}} \\ {{{\Theta^{- 1}\left\lbrack {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I} + {({HA})^{H}{HA}}} \right\rbrack}^{- 1}\begin{matrix} {{({HA})^{H}\mspace{14mu} {for}\mspace{14mu} {unbiased}\mspace{14mu} {Joint}} -} \\ {{MMSE}\mspace{14mu} {Equalizer}} \end{matrix}},} \end{matrix} \right.} & (12) \end{matrix}$

where, Θ¹ is diagonal bias correction matrix for joint-processing, where,

$\Theta = {{diag}{\left\{ {\left\lbrack {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I} + {({HA})^{H}{HA}}} \right\rbrack^{- 1}({HA})^{H}{HA}} \right\}.}}$

Low Complexity GFDM Transmitter:

In this section, low complex GFDM transmitter is presented. A matrix is factorized into special matrices to obtained low complexity transmitter without incurring any assumptions related to GFDM parameters. The stepwise operation of the GFDM transmitter is provided hereunder:

Step 1: Complex valued data symbols corresponding to a particular timeslot and different sub-carriers are modulated to corresponding sub-carrier frequencies using N-point IFFT operation.

Step 2: Modulated data symbols in step 1 are grouped according to sub-carrier numbers, whereby, in each group,

step a: Samples are converted into frequency domain using M-point FFT.

step b: Samples computed in 2(a) are multiplied with a precomputed weight.

step c: Samples in 2(b) are converted back into time domain by using M-point IFFT.

Step 3: Samples obtained after step 2 are regrouped according to time slots.

Signal obtained after step 3 is GFDM transmitted signal.

In the following subsections, the design and implementation of the transmitter is explained.

Low Complexity Transmitter Design

The GFDM modulation matrix A can be given as,

A=P ^(T) U _(M) DU _(M) ^(H) PU _(N),  (13)

where, D=diag{λ(0), λ(1) . . . λ(MN−1)} is diagonal matrix, whose r^(th) element can be given as,

$\begin{matrix} {{\overset{\_}{\lambda}(r)} = {\sum\limits_{m = 0}^{M - 1}\; {{g\left( {{mN} + \left\lfloor \frac{r}{M} \right\rfloor} \right)}{\omega^{m{({r\mspace{11mu} {mod}\mspace{11mu} M})}}.}}}} & (14) \end{matrix}$

and P is a subset of perfect shuffle permutation matrix, which can be defined as, P=[p_(l,q)]0≤l,q≤MN−1, where the matrix element p_(l,q) can be given as,

$\begin{matrix} {P_{l,q} = \left\{ {\begin{matrix} {{1\mspace{14mu} {if}\mspace{14mu} q} = {{lN} + \left\lfloor \frac{1}{M} \right\rfloor}} \\ {0\mspace{14mu} {otherwise}} \end{matrix}.} \right.} & (15) \end{matrix}$

GFDM transmitted signal, x can be given as,

x=P ^(T) U _(M) DU _(M) ^(H) PU _(N) d.  (16)

Lemma 1 Let θ=[θ(0) θ(1) . . . θ(MN−1)]^(T) be a MN length complex valued vector. The vector, {tilde over (θ)}=Pθ=[θ(0) θ(1) . . . {tilde over (θ)}(MN−1)]^(T). The i^(th) element of the vector can be given as,

$\begin{matrix} {{{\overset{\sim}{\vartheta}(i)} = {\vartheta \left( {{\left( {i\mspace{14mu} {mod}\mspace{11mu} M} \right)N} + \left\lfloor \frac{i}{M} \right\rfloor} \right)}},{0 \leq i \leq {{MN} - 1.}}} & (17) \end{matrix}$

The vector, θ=P^(T)θ=[θ(0) θ(1) . . . θ(MN−1)]^(T). The i^(th) element of the vector can be given as,

$\begin{matrix} {{{\overset{\_}{\vartheta}(i)} = {\vartheta \left( {{\left( {i\mspace{14mu} {mod}\mspace{11mu} N} \right)M} + \left\lfloor \frac{i}{N} \right\rfloor} \right)}},{0 \leq i \leq {{MN} - 1.}}} & (18) \end{matrix}$

Low Complexity Transmitter Implementation:

The low complexity transmitter can be obtained using Corollary 1 and Lemma 2. FIG. 3 presents low complexity transmitter implementation.

The vector e=U_(N)d, can be obtained by M, N point IFFT. The vector {tilde over (e)}=Pe can be obtained by shuffling the vector e using (17). The vector c=U^(H) _(M){tilde over (e)} can be obtained using N, M-point FFT's. Using (13), the matrix D, can be pre-computed at the transmitter. The MN length vector, z=Dc, can be obtained by MN-point complex multiplication. The MN length vector, {tilde over (x)}=U_(M)Z can be implemented using N, M-point IFFT. Finally, the transmitted signal, x=P^(T){tilde over (x)} can be obtained by shuffling {tilde over (x)} according to (18).

Low Complexity Two Stage GFDM Receiver:

In this section, the low complexity linear GFDM receivers is disclosed i.e. (1) MF (2) ZF and (3) Biased MMSE and (4) Unbiased MMSE. The stepwise operations of the GFDM receivers are provided hereunder:

step 1: Channel equalized samples are grouped according to sub-carrier numbers. A total N groups are made having M samples and for each group,

step 1a: Samples in a group are converted into frequency domain using M-point FFT;

step 1b: Samples computed in 1a is multiplied with pre-computed weights.

step 1c: Samples in 1(b) are converted back into time domain by using M-point IFFT.

step 2: samples obtained after step 1 are regrouped according to time slots. A total M groups are made having N samples and for each group

step 2a: Samples are converted into frequency domain using N-point FFT;

step 2b: Samples computed in 2a are multiplied with pre-computed weights.

Signal obtained after step 2 is self-interference equalized signal.

Low Complexity GFDM Receiver Design

Receiver in AWGN channel is self-interference equalization. For multipath fading channel, channel equalization is followed by self-interference equalization. Theorem 1 relates to unified low complexity GFDM linear self-interference equalizers. Corollary 1 gives unified implementation of GFDM receivers in AWGN as well as multipath fading channel.

Theorem 1 GFDM equalization matrix A_(eq) can be written in a unified manner as,

A _(eq) =ΘU _(N) ^(H) P ^(T) U _(M) D _(eq) U _(M) ^(H) P,  (19)

where, D_(eq) is a diagonal MN-order matrix, which can be given as,

$D_{eq} = \left\{ \begin{matrix} {D^{H}\mspace{11mu} {for}\mspace{14mu} {MF}\mspace{14mu} {Equalizer}} \\ {D^{- 1}\mspace{11mu} {for}\mspace{14mu} {ZF}\mspace{14mu} {Equalizer}} \\ {\left\lbrack {R_{v} + {D^{H}D}} \right\rbrack^{- 1}D^{H}} \\ {{for}\mspace{14mu} {biased}\mspace{14mu} {MMSE}\mspace{14mu} {Equalizer}} \\ {{\Theta \left\lbrack {R_{v} + {D^{H}D}} \right\rbrack}^{- 1}D^{H}} \\ {{{for}\mspace{14mu} {unbiased}\mspace{14mu} {MMSE}\mspace{14mu} {Equalizer}},} \end{matrix} \right.$

and Θ=Θ_(gfdm) ⁻¹ for unbiased MMSE and Θ=I_(MN) for other equalizers. Further, Θ_(gfdm) can be given as,

$\Theta_{gfdm} = {\frac{1}{MN}{\sum\limits_{r = 0}^{{MN} - 1}\; {\frac{{\lambda_{r}}^{2}}{{\lambda_{r}}^{2} + \frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}}{I_{MN}.}}}}$

Corollary 1 The estimated data, d, can be given as,

$\begin{matrix} {d = \left\{ \begin{matrix} {\Theta \; U_{N}^{H}P^{T}U_{M}D_{eq}U_{M}^{H}{Pz}\mspace{14mu} {for}\mspace{14mu} {AWGN}\mspace{14mu} {Channel}} \\ {\Theta \; U_{N}^{H}P^{T}U_{M}D_{eq}U_{M}^{H}{PW}_{MN}\Lambda_{eq}W_{MN}^{H}z} \\ {{for}\mspace{14mu} {Multipath}\mspace{14mu} {Fading}\mspace{14mu} {Channel}} \end{matrix} \right.} & (20) \end{matrix}$

Low Complexity Receiver Implementation:

The low complexity structure of GFDM self-interference cancellation can be obtained by using Corollary 1 and Lemma 1.

Channel Equalization.

To implement y¹=Λ_(eq)W_(MN) ^(H)z, MN-point FFT of z is multiplied with Λ_(eq). Finally, MN-point IFFT of y¹ is taken to implement y=W_(MN)y¹.

Self-interference Equalization.

The vector y=Py, can be obtained by shuffling the y vector using (16). The MN×1 vector α=U_(M) ^(H)y can be implemented by using N, M-point IFFT's. The vector α is then multiplied to the diagonal matrix D_(eq) to obtain β. The vector θ=U_(M)β can be implemented using N, M-point FFTs. The vector, {tilde over (θ)}=P^(T)θ, can be implemented by shuffling the θ vector using (17). Now, the vector, d=U_(N){tilde over (θ)} can be implemented using M, N-point FFTs. Finally, d=Θd can be obtained by using MN-point multiplier.

Low complexity Joint-MMSE Receiver:

Step 1: Received signal is converted into frequency domain using M N-point FFT.

Step 2: Received signal in frequency domain is multiplied with complex valued channel information-based weights and then converted back to time domain using MN-point IFFT.

Step 3: Samples obtained in step 2 are grouped according to subcarrier number. N such groups are formed having M samples each and for each group

Step 3a: Samples are converted into frequency domain using M-point FFT and multiplied with pre-computed weights;

Step 4 Samples obtained in step 3a are regrouped according to time slot number. M such groups are formed having N samples and for each group

Step 4a: Samples are converted into frequency domain using N-point FFT.

Step 4b: Samples obtained in step 4a are processed as explained in Algorithm 1 or Algorithm 2.

Step 5 Samples obtained after step 4 are regrouped according to sub-carrier number. N such groups are formed having M samples each and for each group, repeat

Step 5a Samples are converted back to time domain using M-point IFFT and multiplied with complex weights which are computed using Algorithm 3.

Step 6 Samples obtained after step 5 are regrouped according to time slots to obtain equalized MN-point samples.

Joint MMSE Low Complexity Receiver Design

Theorem 1 The estimated data vector for joint-MMSE receiver, d, can be given as,

$\begin{matrix} \begin{matrix} {d = {\Theta_{J\; P}^{- 1}F_{b}{EU}_{N}^{H}D^{H}F_{b}^{H}W_{MN}\Lambda^{H}W_{MN}^{H}y}} \\ {= {\Theta_{J\; P}^{- 1}P^{T}U_{M}{PEU}_{N}^{H}D^{H}P^{T}U_{M}^{H}{PW}_{MN}\Lambda^{H}W_{MN}^{H}{y\left( {{UsingLemma}\; 1} \right)}}} \end{matrix} & (21) \end{matrix}$

where, E=diag{E₀, E₁, . . . E_(M-1)}, is a MN×MN size block diagonal matrix with blocks of size N×N. A block of E,

$E_{u} = \left\lbrack {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I} + C_{u}} \right\rbrack^{- 1}$

where, 0≤u≤M−1, of size N×N, where, Θ_(JP)=Θ for unbiased receiver and I_(MN) for biased MMSE receiver. Further, C_(u) can be given as,

C _(u) =L _(u) ^(H)γ_(u) L _(u), 0≤u≤M−1,  (22)

where, L_(u) is a circulant N×N matrix which can be represented in terms of its first column as, L_(u)=circ{l_(u)(0), l_(u)(1), . . . l_(u)(N−1)}, where, p^(th) element of the first column can be given as,

${{l_{u}(p)} = {\frac{1}{N}{\sum\limits_{r = 0}^{N - 1}\; {{\lambda \left( {{uN} + r} \right)}e^{\frac{j\; 2\; {\pi {({- p})}}r}{N}}e^{\frac{{- j}\; 2\; \pi \; {ur}}{MN}}}}}},$

0≤u≤M−1, and, γ_(u) is a N×N diagonal matrix which can be given as,

γ_(u)=diag{{tilde over (h)}(u)|² ,|{tilde over (h)}(M+u)|² . . . |{tilde over (h)}((N−1)M+u)|²},0≤u≤M−1.  (23)

Further, Θ can be given as,

Θ=I _(M) ⊗S,  (24)

where,

${S = {\frac{1}{M}{\sum\limits_{u = 0}^{M - 1}\; {{diag}\left\{ {\overset{\sim}{E}}_{u} \right\}}}}},$

where,

${\overset{\sim}{E}}_{u} = {\left\lbrack {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I_{N}} + C_{u}} \right\rbrack^{- 1}C_{u}}$

is a N×N matrix, 0≤u≤M−1. Proof. The theorem can be proved using the fact that (H A)^(H) H A is block circulant matrix with blocks of size N×N.

The implementation of joint-MMSE receiver based on (15), requires the inversion of M number of, N×N size

$\left\lbrack {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I} + C_{u}} \right\rbrack$

matrices. The direct implementation of this inversion requires O(MN³). To see the possibility of further reduction in complexity, E_(u) is factorized for 0≤u≤M−1, Factorization of E_(u) for 0≤u≤M−1:

Using (16), E_(u) can be written as,

${E_{u} - \left\lbrack {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I} + {L_{u}^{H}\mathrm{\Upsilon}_{u}L_{u}}} \right\rbrack^{- 1}} = {{L_{u}^{- 1}\left\lbrack {\mathrm{\Upsilon}_{u} + {\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}\left( {L_{u}L_{u}^{H}} \right)^{- 1}}} \right\rbrack}^{- 1}{\left( L_{u}^{H} \right)^{- 1}.}}$

Since L_(u) is a circulant matrix, it can be further factorized as,

L _(u) =W _(N) R _(u) W _(N) ^(H),

where, R_(u) is a diagonal matrix of order N, which can be computed as, R_(u)=W_(N) ^(H)L_(u)W_(N). Using this, (19) can be written as,

E _(u) =W _(N) R _(u) ⁻¹ W _(N) ^(H)Φ_(u) W _(N)(R _(u) ^(H))⁻¹ W _(N) ^(H),

where,

$\Phi_{u} = {℧_{u}^{- 1} = {\left\lbrack {\mathrm{\Upsilon}_{u} + {\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}\left( {L_{u}L_{u}^{H}} \right)^{- 1}}} \right\rbrack^{- 1}.}}$

Using the properties of circulant matrices, it can be shown that L_(u)L_(u) ^(H) is a circulant matrix because L_(u) is a circulant matrix. Further, L_(u)L_(u) ^(H)=W_(N)|R_(u)|²W_(N) ^(H). Using this, (L_(u)L_(u) ^(H))⁻¹=W_(N)|R_(u)|⁻²W_(N) ^(H). Hence,

_(u) is a circulant-plus-diagonal matrix. It can be easily seen that elements of γ_(u) and |R_(u)|⁻² are positive. It can be concluded that

_(u) is a positive definite matrix too.

Since

_(u) is a positive definite matrix, inversion of

_(u) can be computed using Conjugate Gradient (CG) algorithm. CG algorithm gives exact solution in N iterations. Hence a direct implementation of joint-MMSE receiver can be obtained using CG method. In each iteration, a matrix-vector multiplication is required. In our case this matrix is a circulant-plus-diagonal matrix. Using the properties of circulant matrix, matrix-vector multiplication can be implemented using N-point FFT and IFFT. Thus, direct implementation of joint-MMSE receiver requires O(MN² log₂N) computations.

It has been showed that to implement the receiver, the most computationally complex operation is to invert N-order

_(u) matrix, 0≤u≤M−1. First the structure of

_(u) matrix is investigated. The low complexity bias correction is also investigated.

Structure of

_(u) matrix, 0≤u≤M−1:

Using the properties of circulant matrices, it can be shown that diagonal values of (L_(u)L_(u) ^(H))⁻¹ are equal and can be given as, diag{(L_(u)L_(u) ^(H))⁻¹}=μ_(u)I_(N), where, μ_(u) can be given as,

$\mu_{u} = {\frac{\sum\limits_{r = 0}^{N - 1}\; \left\lbrack {R_{u}}^{- 2} \right\rbrack_{r,r}}{N}.}$

Next, two matrices are defined,

$\Delta_{u} = {\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}\left\{ {\left( {L_{u}L_{u}^{H}} \right)^{- 1} - {\mu_{u}I_{N}}} \right\}}$ and $_{u} = {\mathrm{\Upsilon}_{u} + {\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}\mu_{u}{I_{N}.}}}$

Using this,

_(u) can be written as,

_(u)=[

_(u)+Δ_(u)].

It can be seen that

_(u) is a diagonal matrix and Δ_(u) is a circulant matrix with zero diagonal values. Using the properties of circulant matrix, ∥Δ_(u)∥² can be given as,

${\Delta_{u}}^{2} = {\left( \frac{\sigma_{v}^{2}}{\sigma_{d}^{2}} \right)^{2}{trace}{\left\{ R_{u}^{2} \right\}.}}$

E[∥

_(u)∥²] can be approximated as,

${E\left\lbrack {_{u}}^{2} \right\rbrack} \approx {{N\left( {\left( {1 + {\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}\mu_{u}}} \right)^{2} + 1} \right)}.}$

It is observed that E[∥

_(u)∥²] is approximately independent of assumed channel power delay profile.

Next, the ratio of power in Δ_(u) to the power in

_(u) is analyzed. To do so, ρ=

$\frac{1}{M}{\sum\limits_{u = 0}^{M - 1}\; {E_{h}\left\lbrack \frac{\Delta_{n}}{_{u}} \right\rbrack}}$

is defined. Urban micro (UMi) channel is taken and Monte-Carlo simulations are persomed to compute

${E_{h}\left\lbrack \frac{\Delta_{u}}{_{u}} \right\rbrack}.$

It is further averaged over u to obtain ρ. Plot of ρ versus SNR is given in FIG. 5A for N∈{32,128,512}, M=3 and ROF∈{(0.1,0.9}. It is observed that ρ decreases with SNR and increases with ROF. For same SNR and ROF value, it decreases with increasing value of N. FIG. 5B plots ρ against SNR for M∈{3,11,21}, N=16 and ROF∈{(0.1,0.9}. It can be seen that ρ decreases with SNR and increases with ROF value. For low ROF values, ρ is decreases with M and for high ROF values, it increases with M. To see the effect of h, three different channel lengths of following uniform power delay profile (PDP) is considered. The value of ρ is plotted against SNR in FIG. 5C for channel lengths L=[16,32,64] representing different rms delay spread σ_(rms)≈[1,1.8,3.7] μs respectively. It can be observed that ρ is same for different channel lengths. Hence, it can be concluded that power in

_(u) is much higher than the power in Δ_(u) irrespective of channel propagation conditions.

$\psi = {E_{h}\left\lbrack \frac{\Delta_{n}}{_{u}} \right\rbrack}$

is defined. In FIG. 5D, the value of ψ against u∈[0M−1] and ROF=0.3,0.5,0.9 for N=64, M=5,15 and SNR=10 dB is plotted. It is observed that value of ψ is highest for middle two values of u. In case of M=15 and ROF=0.9, the power in Δ_(u) for u=7,8 is comparable to the power in

_(u). Taylor Series Method for Computing Φ_(u)=

_(u), 0≤u≤M−1:

It can be concluded from the previous section that the power in

_(u) is much higher than the power in Δ_(u) when both M and ROF have value. Therefore, Taylor series expansion of Φ_(u) can be used for other situations. In such cases, Φ_(u)≈

_(u) ⁻¹−

_(u)Δ_(u)

_(u) ⁻¹+

_(u) ⁻¹Δ_(u)

_(u) ⁻¹Δ_(u)

_(u) ⁻¹+ . . . . Also, using (19) and (21), the following expression can be reached

Δ_(u)=W_(N)R_(u)W_(N) ^(H), where

$R_{u} = {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}{R_{u}}^{- 2}} - {\mu_{i}{I_{N}.}}}$

Therefore,

Φ_(u)≈

_(u) ⁻¹[I _(N) −W _(N) R _(u) W _(N) ^(H)

_(u) ⁻¹ +W _(N) R _(u) W _(N) ^(H)

_(u) ⁻¹ W _(N) R _(u) W _(N) ^(H)

_(u) ⁻¹+ . . . ].  (25)

Algorithm to Multiply E in (14) with Complex Valued Vector Using Taylor Series Method:

κ=U_(N) ^(H)D^(H)F_(b) ^(H)W_(MN)Λ^(H)W_(MN) ^(H)y, is considered to be an intermediate MN length complex valued vector in (14). Now to compute ρ=Eκ, another intermediate vector in (14), E can computed by putting (23) in (20). Algorithm 1 discuss a low complexity algorithm to multiply ρ=Eκ using Taylor Series expansion. This algorithm is used to develop low complexity receiver structure. Matrices vectors and scalars which are only related GFDM parameters (do not depend on the channel) can be precomputed at the receiver. Hence, it is assumed that, the knowledge of R_(u), |R_(u)|⁻², (R_(u) ^(H))⁻¹, R_(u) ⁻¹, R_(u), μ_(u) for 0≤u≤M−1. Also, it is assumed that the knowledge of

$\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}$

which can be computed at the receiver and which is known from the channel estimation. K_(T)+1 is the number of terms in (24) i.e. the number of iterations for step 7 to step 9 is K_(T).

Algorithm 1 Computation of 

 = Eκ Using Taylor Series Method    1: ${{Precomputed}:R_{u}},{R_{u}}^{- 2},\left( R_{u}^{H} \right)^{- 1},R_{u}^{- 1},{\overset{\sim}{R}}_{u},µ_{u},{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}\mspace{14mu} {and}}$ Λ for 0 ≤ u ≤ M − 1  2: for u = 0 : M − 1 do  3:  Compute : Υ_(u) using (16)  4:   ${{Compute}:_{u}} = {\mathrm{\Upsilon}_{u} + {\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}µ_{u}I_{N}}}$  5:  Compute : 

 = W_(N)(R_(u) ^(H))⁻¹W_(N) ^(H)κ_(u)  6:  Take : t⁽⁰⁾ = 

   7:  for k = 1 : K_(T) do  8:   t^((k)) = W_(N){tilde over (R)}_(u)W_(N) ^(H) 

_(u) ⁻¹t^((k−1))  9:  end for 10:   ${{Compute}:\varrho_{u}^{(1)}} = {_{u}^{- 1}{\sum\limits_{k = 0}^{K_{T}}t^{(k)}}}$ 11:  Compute : 

_(u) = W_(N)R_(u) ⁻¹W_(N) ^(H) 

_(u) ⁽¹⁾ 12: end for 13:

 = [ 

₀ 

₁ . . . 

_(M−1)]^(T) 14: return 

  Low Iteration Conjugate Gradient (CG) Method for Computing Φ_(u)=

_(u) ⁻¹, 0≤u≤M−1

It is established earlier, that

_(u) is a positive definite matrix. So,

_(u)

_(u)=κ_(u) can be computed using CG algorithm. CG algorithm gives exact solution in N iterations. To reduce the complexity further Jacobi precoded CG method is used i.e. the system

_(u) ⁻¹

_(u)=

_(u) ⁻¹κ_(U) is solved. Thus, [I_(N)+

Δ_(u)]

_(u)=

_(u) ⁻¹κ_(u) is solved using CG method. Different iteration count is also used for different values of u since

$\frac{\Delta_{u}}{_{u}}$

changes with u (as illustrated in FIG. 5D). When

$\frac{\Delta_{u}}{_{u}}$

is small, number of iterations can be kept small. When

$\frac{\Delta_{u}}{_{u}}$

is large, number of iteration is made large to obtain low errors. Algorithm to Multiply E in (14) with Complex Valued Vector Using CG Method:

Algorithm 2 discuss a low complexity algorithm to multiply ρ=Eκ using CG method. This algorithm is used to develop low complexity receiver structure. k_(C) is considered to be a M-length vector which holds iteration counts for different values of u. Same as in Algorithm 1, matrices, vectors and scalars which are only related GFDM parameters (do not depend on the channel) are assumed to be precomputed at the receiver.

Algorithm 2 Computation of 

 = Eκ using conjugate gradient (CG) method    1: ${{Precomputed}:R_{u}},{R_{u}}^{- 2},\left( R_{u}^{H} \right)^{- 1},R_{u}^{- 1},{\overset{\sim}{R}}_{u},µ_{u},{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}\mspace{14mu} {and}}$ Λ for 0 ≤ u ≤ M − 1  2: for u = 0 : M − 1 do  3:  Compute : Υ_(u) using (17)  4:   ${{Compute}:_{u}} = {\mathrm{\Upsilon}_{u} + {\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}µ_{u}I_{N}}}$  5:  Compute : 

 = W_(N)(R_(u) ^(H))⁻¹W_(N) ^(H)κ_(u)  6:  Compute : 

⁽¹⁾ = 

_(u) ⁻¹ 

 7:  Take : w(0) = 

⁽¹⁾  8:  Compute : e(0) = 

⁽¹⁾ − [I_(N) + 

_(u) ⁻¹W_(N){tilde over (R)}_(u)W_(N) ^(H)]w(0)  9:  Take : p(0) = e(0) 10:  for k = 0 : k_(C)(u) do 11:   Compute : α(k) = [I_(N) + 

_(u) ⁻¹W_(N){tilde over (R)}_(u)W_(N) ^(H)]p(k) 12:    ${{Compute}:{℧(k)}} = \frac{{e^{H}(k)}{e(k)}}{{p^{H}(k)}{\alpha (k)}}$ 13:   Compute : w(k + 1) = w(k) + 

(k)p(k) 14:   Compute : e(k + 1) = e(k) − 

(k)α(k) 15:    ${{Compute}:{\tau (k)}} = \frac{{e^{H}\left( {k\text{+}1} \right)}{e\left( {k\text{+}1} \right)}}{{e^{H}(k)}{e(k)}}$ 16:   p(k + 1) = e(k + 1) + τ(k)p(k) 17:  end for 18:  Take : 

_(u) ⁽¹⁾ = w(k_(C)(u)) 19:  Compute : 

_(u) = W_(N)R_(u) ⁻¹W_(N) ^(H) 

_(u) ⁽¹⁾ 20: end for 21:

 = [ 

₀ 

₁ . . . 

_(M−1)]^(T) 22: return 

 

Approximation of Bias Correction Matrix Θ

To compute Θ, diag{{tilde over (E)}_(u)}, 0≤u≤M−1, is required (see (17)). Let, Q_(u)=diag{C_(u)} and S_(u)=C_(u)−Q_(u), using (17) {tilde over (E)}_(u) can be given as,

${\overset{\sim}{E}}_{u} = {{\left\lbrack {\underset{J_{u}}{\underset{}{{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I_{N}} + Q_{u}}} + S_{u}} \right\rbrack^{- 1}\left\lbrack {Q_{u} + S_{u}} \right\rbrack} = {{\left\lbrack {J_{u} + S_{u}} \right\rbrack^{- 1}\left\lbrack {Q_{u} + S_{u}} \right\rbrack}.}}$

It can be shown that ∥Q∥>>∥S_(u)∥. This implies that ∥J_(u)∥>>∥S_(u)∥. Hence, [J_(u)+S_(u)]⁻¹ can be approximated using Taylor series. {tilde over (E)}_(u) can be approximated as, {tilde over (E)}_(u)≈[J_(u) ⁻¹−J_(u) ⁻¹S_(u)J_(u) ⁻¹][Q_(u)+S_(u)]. Using this diag{{tilde over (E)}_(u)} can be approximated as, diag{{tilde over (E)}_(u)}≈J_(u) ⁻¹Q_(u)−diag{J_(u) ⁻¹S_(u)J_(u) ⁻¹S_(u)}. Since ∥J_(u)∥>>∥S_(u)∥ as well as ∥Q_(u)∥>>∥S_(u)∥, it can be easily shown that ∥J_(u) ⁻¹Q_(u)∥>>∥diag{J_(u) ⁻¹S_(u)J_(u) ⁻¹S_(u)}∥. So, diag{{tilde over (E)}_(u)} can be further approximated as,

diag{{tilde over (E)} _(u) }≈J _(u) ⁻¹ Q _(u)

Now, Q_(u)=diag{C_(u)} to be computed for computation of diag{{tilde over (E)}_(u)} in (25). It can be shown, that only three elements in a column of L_(u) matrix is dominant, which are l_(u)(0), l_(u)(1) and l_(u)(N−1). Other elements have comparatively lesser power by at least 40 dB. Hence L_(u) matrix is approximated as, L_(u)=circ{l_(u)(0), l_(u)(1), 0 . . . 0l_(u)(N−1)}. Using this approximation, (15) and (16), the diagonal elements of Q_(u) can be approximated as,

$\begin{matrix} {{{{\left\lbrack Q_{u} \right\rbrack_{s,s} \approx {\sum\limits_{{r = 0},1,{N - 1}}^{\;}\; {{l_{u}(r)}}^{2}}}}}{\overset{\sim}{h}\left( {u + {\left( {\left( {r + s} \right){mod}\; N} \right)M}} \right)}{^{2},}} & (26) \end{matrix}$

where 0≤s≤N−1 and 0≤u≤M−1. Now, using (17), (25) and (26), Θ can be approximated as,

$\begin{matrix} {\Theta \approx {\frac{1}{M}{I_{M} \otimes {\sum\limits_{u = 0}^{M - 1}\; {J_{u}^{- 1}{Q_{u}.}}}}}} & (27) \end{matrix}$

Algorithm to Compute Θ_(JP)

Algorithm 3 explains low complexity computation of Θ_(JP). |l_(u)(r)|² is assumed for r=0,1, N−1 and 0≤u≤M−1 is precomputed at the receiver which requires storage of 3M real values for unbiased receiver.

Algorithm 3 Computation of Θ_(JP)    1: if biased receiver then  2:  Θ_(JP) = I_(MN)  3: else  4:   $\quad\begin{matrix} {{{{Precomputed}\mspace{14mu} {{l_{u}(r)}}^{2}\mspace{14mu} {for}\mspace{14mu} r} = 0},1,{N - {1\mspace{14mu} {and}}}} \\ {0 \leq u \leq {M - {1.\mspace{14mu} \frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}\mspace{14mu} {and}\mspace{14mu} \Lambda}}} \end{matrix}$  5:  for u = 0 : M − 1 do  6:   Compute : Q_(u) using (26)  7:    ${{Compute}:J_{u}} = \left\lbrack {{\frac{\sigma_{v}^{2}}{\sigma_{d}^{2}}I_{N}} + Q_{u}} \right\rbrack$  8:   Compute : J_(u) ⁻¹Q_(u)  9:  end for 10:  Compute : Θ_(JP) = Θ using (27) 11: end if 12: return Θ_(JP)

Joint-MMSE Low Complexity Structure

To implement,

y=W _(MN)Λ^(H) W _(MN) ^(H) y

IFFT_(MN) of the product of diagonal matrix Λ^(H) with (FFT_(MN) of y) is computed. To implement, ξ=P^(T)U_(M) ^(H)Py, the vector, y is first shuffled according to (17) and then passed through N FFT_(M) whose output is again shuffled according to (18). To implement, κ=U_(N) ^(H)D^(H)ξ, the vector, ξ, is multiplied with D^(H) using MN-point multiplier whose output is then passed through M FFT_(N). The vector, κ is then passed through M N-order square matrix inversion block to obtain ρ=Eκ using Algorithm 1 or 2. In the last, the vector, κ is first shuffled according to (17) and then passed through N, IFFT_(M) whose output is again shuffled according to (18) and multiplied to Θ_(JP) to obtain estimated data, d. Θ_(JP) can be computed using Algorithm 3.

Testing Two Stage Receiver

Monte-Carlo simulation is performed for GFDM system which comprises of the proposed transmitter and two-stage receiver. Each point in the BER curve is calculated for 107 transmission bits.

BER of the proposed low complexity transceiver in multipath fading channel is plotted in FIG. 7. The proposed receivers do not incur any performance loss over direct implementations. The proposed MMSE receiver shows significant BER performance gain over other GFDM receivers. For ROF value of 0.1, the proposed MMSE receiver gives the best performance. MMSE receiver achieves SNR gain of 4 dB over OFDM at BER of 10-5. This SNR gain over OFDM is due to higher frequency resolution of GFDM [Ref: N. Michailow, M. Matthe, I. Gaspar, A. Caldevilla, L. Mendes, A. Festag, and G. Fettweis, \Generalized Frequency Division Multiplexing for 5th Generation Cellular Networks,” IEEE Transactions on Communications, vol. 62, no. 9, pp. 3045{3061, September 2014]. MMSE receiver also achieves SNR gain of 2.5 dB over our proposed ZF receiver and SIC receiver in [Ref: A. Farhang, N. Marchetti, and L. E. Doyle, \Low complexity GFDM receiver design: A new approach,” in 2015 IEEE International Conference on Communications (ICC), June 2015, pp. 4775-4780.] at BER of 4×10-6. For ROF value of 0.9, the proposed MMSE receiver shows a SNR gain of 5 dB over our proposed ZF receiver at BER of 10-4. BER of SIC receiver in [Ref: Ref: A. Farhang, N. Marchetti, and L. E. Doyle, \Low complexity GFDM receiver design: A new approach,” in 2015 IEEE International Conference on Communications (ICC), June 2015, pp. 4775-4780.] floors at 10-2 and has SNR loss of 15 dB over the proposed MMSE receiver.

TABLE 1 Simulation Parameter for two stage receiver Parameters Value Number of Subcarriers N 128 Number of Timeslots M  8 Mapping 16 QAM Pulse shape Raise Cosine (RC) Sub-carrier Bandwidth 15 KHz Channel Extended Typical Urban (ETU) [Ref: M. Series, “Guidelines for evaluation of radio interface technologies for IMT- advanced”, Report ITU, no. 2135-1, 2009.] Carrier Frequency 2.4 GHz Maximum Doppler shift 100 Hz RMS delay Spread 1 micro second Channel Equalization MMSE FDE

Joint MMSE Receiver

In this section, BER performance of proposed receiver is presented in the multipath channel. Table 8 presents the simulation parameters. The multipath channel Urban Micro (UMi) [26] with 20 taps is considered, whose channel delay and channel power are [0 10 15 20 35 40 45 50 55 200 205 250 330 440 515 530 580 590 625 730] ns and [−6.7 −4.9 −7.1 −1.9 −6.3 −3 −5.2 −7 −7.5 −10.8 −5.2 −4.9 −9.2 −15.5 −12.4 −16.9 −12.7 −23.5 −22.1 −23.6] dB, respectively. The CP is chosen long enough to accommodate the wireless channel delay spread. A coded system with code rate of 0.5 is assumed. Convolution code is used with constraint length of 7 and code generator polynomial of 171 and 133. A random interleaver having length equal to K_(qam)MN is considered, where K_(qam) is number of bits in a QAM symbol. Soft maximum likelihood (ML) decoding is implemented at the receiver. Each point in our BER curve is calculated for 10⁸ transmission bits.

TABLE 2 Simulation Parameters for Joint-MMSE Receiver Parameters Value Number of Subcarriers N 128 Number of Timeslots M  5 Mapping 16 QAM Pulse shape Raise Cosine (RC) Sub-carrier Bandwidth 15 KHz Channel Umi [Ref: M. Series, “Guidelines for evaluation of radio interface technologies for IMT-advanced”, Report ITU, no. 2135-1, 2009.] Carrier Frequency 2.4 GHz Maximum Doppler shift 6.67 Hz RMS delay Spread 0.38 micro second Decoding Soft Maximum Likelihood (ML)

For soft ML decoding, post processing SNR of MMSE receiver output is required. It is straight forward to compute SNR (Γ(1)) for l^(th) symbol which can be given as,

${\Gamma (l)} = {\frac{1}{\lbrack\Theta\rbrack_{l,l} - 1}.}$

As discussed, Θ is computed for correcting the bias of MMSE equalization outputs using Algorithm 3. Θ is periodic with N. Thus Γ(l+mN)=Γ(l), where, m∈[0, M−1]. This means that computation of F requires additional N complex multiplications.

The BER performance of the proposed receiver is computed with the direct ones in Michailow et al., “Generalized Frequency Division Multiplexing for 5th Generation Cellular Networks.” for N=128, M=5 i.e. in TI scenario for ROF value of 0.3 and 0.9 in FIGS. 8A and 8B respectively. For CG based method, k_(C)=[1410104] is considered i.e. K_(C)=5. When K_(T)=0, the proposed TS based receiver shows 0.5 dB SNR loss at the BER of 2×10⁻⁶ for ROF value of 0.3. Whereas for ROF value of 0.9, the proposed receiver with, K_(T)=0, has 1 dB SNR degradation at the BER of 10⁻⁵. When K_(T)=2 and K_(C)=5, It is observed that there is no notable degradation in BER as compared to direct implementation for ROF values of 0.3 as well as 0.9. The bias correction provides SNR gain of 1 dB and 0.8 dB at BER of 10⁻⁵ for ROF value of 0.3 and 0.9 respectively without any observable increase in complexity. It can be concluded that approximations made in Algorithm 1 and Algorithm 2 have no noticeable effect on BER performance. The proposed bias correction has BER benefits without any observable increase in complexity. 

What is claimed is:
 1. A generalized frequency division multiplexing (GFDM) transceiver system comprising: low complex GFDM transmitter with multiple sub-carriers and timeslots having IFFT based modulator for modulating data corresponding to a particular timeslot and different sub-carriers to corresponding sub-carrier frequencies and thereby generating transmittable GFDM data signal; multipath frequency selective fading GFDM channel having uncorrelated channel coefficients corresponding to different paths for transmitting the modulated GFDM data signal; and low complex GFDM receiver configured to operate with said multipath frequency selective fading channel involving channel equalization followed by self-interference equalization to receive the transmitted modulated GFDM data signal and thereby de-modulate the GFDM data signal to obtain the data.
 2. The generalized frequency division multiplexing (GFDM) transceiver system as claimed in claim 1, wherein the low complex GFDM transmitter includes: N-point IFFT operator to receive the data corresponding to a particular timeslot and different sub-carriers and modulate the same to corresponding sub-carrier frequencies; means for shuffling physical connections in the N-point IFFT operator for grouping the modulated data to sub-carrier numbers, whereby, in each group, the data is converted into frequency domain using M-point FFT operator and multiplied with a precomputed weight and thereafter converted back into time domain by using M-point IFFT operator; and means for shuffling physical connections in the M-point IFFT operator for grouping the data according to time slots and generate transmittable GFDM data signal.
 3. The generalized frequency division multiplexing (GFDM) transceiver system as claimed in claim 1, wherein the low complex GFDM transmitter generated transmittable GFDM data signal is ${x\lbrack n\rbrack} = {\frac{1}{\sqrt{N}}{\sum\limits_{m = 0}^{M - 1}\; {\sum\limits_{k = 0}^{N - 1}\; {d_{m,k}{g\left\lbrack {n - {mN}} \right\rbrack}_{MN}e^{\frac{j\; 2\pi \; {nk}}{N}}}}}}$ for N sub-carriers and M timeslots, where g(n), n=0, 1, . . . , MN−1 is MN length filter response and d_(m,k)∈C, m=0, 1, . . . , M−1, k=0, 1, . . . , N−1 is QAM modulated data symbol.
 4. The generalized frequency division multiplexing (GFDM) transceiver system as claimed in claim 3, wherein the GFDM transmitted signal is critically sampled Inverse Discrete Gabor Transform (IDGT) of d by using the IDGT matrix factorization whereby Modulation Matrix, A can be given as: $A = {{G \times U_{N}} = {\begin{pmatrix} \Psi_{0} & \Psi_{M - 1} & \ldots & \Psi_{1} \\ \Psi_{1} & \Psi_{0} & \ldots & \Psi_{2} \\ \vdots & \vdots & \ddots & \vdots \\ \Psi_{M - 1} & \Psi_{M - 2} & \ldots & \Psi_{0} \end{pmatrix}\begin{pmatrix} W_{N} & \; & \; \\ \; & \ddots & \; \\ \; & \; & W_{N} \end{pmatrix}}}$ where, Ψ_(m)=diag{g[mN], g[mN+1], . . . , g[mN+N−1]} for 0≤m≤M−1, is N×N diagonal matrix and W_(N) is N×N normalized IDFT matrix.
 5. The generalized frequency division multiplexing (GFDM) transceiver system as claimed in claim 1, wherein the uncorrelated channel coefficients corresponding to different paths for transmitting the modulated GFDM data signal constitutes channel impulse response vector given as h=[h₀, h₁, . . . h_(L-1)]^(T) where L is channel length and h_(i), for 0≤i≤L−1, represents complex baseband channel coefficient of (i+1)^(th) path, which is assumed to be zero mean circular symmetric complex Gaussian whereby received vector of length N_(CP)+NM+L−1 (for N_(cp)≥L) is given by, Z _(cp) =h*X _(cp) +v _(cp) where v_(cp) is AWGN vector of length MN+N_(CP)+L−1 with elemental variance 6²v.
 6. The generalized frequency division multiplexing (GFDM) transceiver system as claimed in claim 1, wherein the low complex GFDM receiver includes two staged receiver or Joint-MMSE Receiver whereby, the data obtained from the received GFDM data signal by involving channel equalization followed by self-interference equalization.
 7. The generalized frequency division multiplexing (GFDM) transceiver system as claimed in claim 6, wherein the two-staged receiver includes two staged receiver includes: M-point FFT operator for grouping channel equalized received GFDM data signal according to sub-carrier numbers followed by: converting samples of each group into frequency domain using the M-point FFT operator; multiplying the converted samples with pre-computed weights by multiplier means; and converting back the multiplied samples into time domain by using M-point IFFT operator; means for shuffling physical connections in the M-point IFFT operator for regrouping the converted samples according to time slots followed by converting samples of each group into frequency domain using N-point FFT operator; and multiplying the converted samples with pre-computed weights by multiplier means to obtain self-interference equalized data signal.
 8. The generalized frequency division multiplexing (GFDM) transceiver system as claimed in claim 6, wherein the joint-MMSE receiver includes: MN-point FFT operator to convert the equalized received GFDM data signal into frequency domain; means to multiply the received signal in frequency domain with complex valued channel information-based weights and then converting back to time domain using MN-point IFFT operator; means for reshuffling the physical connections in the MN-point IFFT operator for grouping the time domain converted samples according to subcarrier number followed by: converting the samples of each group into frequency domain using the M-point FFT operator; and multiplying the converted samples with pre-computed weights; means for reshuffling the physical connections in the M-point FFT operator for regrouping the converted samples according to time slots followed by converting samples of each group into frequency domain using the N-point FFT operator; processing the converted samples following Algorithm 1; means for reshuffling the physical connections in the N-point FFT operator for regrouping the processed samples according to sub-carrier number followed by: converting the samples of each group to time domain using the M-point IFFT operator; and multiplying the converted samples with complex weights computed using Algorithm 2; and means for reshuffling the physical connections in the M-point IFFT operator for regrouping the multiplied samples according to time slots to obtain equalized MN-point samples.
 9. The generalized frequency division multiplexing (GFDM) transceiver system as claimed in claim 8, wherein the Algorithm 1 enabling development of the low complexity receiver structure including low complexity multiplication to obtain ρ=Eκ using Taylor Series expansion.
 10. The generalized frequency division multiplexing (GFDM) transceiver system as claimed in claim 9, wherein the Algorithm 2 enabling development of the low complexity receiver structure including low complexity multiplication to obtain ρ=Eκ using CG method. 